UG-FT-247/09 
CAFPE-117/09 



Cosmic ray knee and new physics at the TeV scale 

Roberto Barcelo^, Manuel Masip^, lacopo Mastromatteo^ 

^ CAFPE and Departamento de Fisica Tedrica y del Cosmos 
Universidad de Granada, E- 18071 Granada, Spain 

"^International School for Advanced Studies (SISSA) 
Via Beirut 2-4, 1-34014 Trieste, Italy 

rbarcelo@ugr.es, masip@ugr.es, mastroma@sissa.it 



Abstract 

We analyze the possibility that the cosmic ray knee appears at an energy threshold 
where the proton-dark matter cross section becomes large due to new TeV physics. It 
has been shown that such interactions could break the proton and produce a diffuse 
gamma ray flux consistent with MILAGRO observations. We argue that this hypothesis 
implies knees that scale with the atomic mass for the different nuclei, as KASKADE 
data seem to indicate. We find that to explain the change in the spectral index in 
the flux from E^'^-'^ to E^^'"^ the cross section must grow like £'0-^+^ above the knee, 
where j3 = 0.3-0.6 parametrizes the energy dependence of the age (r oc E~l^) of the 
cosmic rays reaching the Earth. The hypothesis also requires mbarn cross sections 
(that could be modelled with TeV gravity) and large densities of dark matter (that 
could be clumped around the sources of cosmic rays). We argue that neutrinos would 
also exhibit a threshold sX E = (m^/mp) Ef^^ee ~ 10^ GeV where their interaction 
with a nucleon becomes strong. Therefore, the observation at ICECUBE or ANITA of 
standard neutrino events above this threshold would disprove the scenario. 



1 Introduction 



The observed cosmic ray flux reaching the Earth extends up to energies around 10^^ GeV. 
It seems remarkable that over 10 decades of energy this flux can be described in very simple 
terms [1]. Between 10 and 10^ GeV (the knee) it is given by 

^«1.8g-" 7'—^ . (1) 

ah/ cm"^ s sr GeV 

At Eknee the spectral index changes to —3.1 and stays constant up to ~ 10^'^ GeV (the 
ankle). There it goes back to —2.7 until (arguably) hits the GZK suppression a decade later 
(see [2], [3] for a review). 

The ankle is generally explained as the overtaking in the cosmic ray flux of a new com- 
ponent of extragalactic origin. These ultra-high energy particles arrive from all regions in 
the sky, which excludes the possibility that they are produced only in the galactic disc. For 
less energetic particles, E < 10^ GeV, one can assume production in the disc despite the fact 
that the observed flux is also isotropic. The reason is that their Larmor radius is smaller, so 
the trajectory is isotropized by the random magnetic fields present in our galaxy [r^ ^ 0.1 
kpc for a. E = 10^ GeV proton inside a. B = 1 fiG magnetic field). 

The other main feature in the spectrum, the knee, seems more involved. It is usually 
thought that E^nee reflects the maximum energy reached by the dominant cosmic accelerators 
in our galaxy. In that case, however, one would expect a transition regime (not seen) between 
Eknee and the energy E where the new acceleration mechanism dominates. Propagation and 
confinement within the galaxy could also play an important role. In particular, the knee 
could correspond to a critical energy where a large fraction of cosmic rays escapes the galactic 
disc before reaching the Earth [1]. Again, it seems difficult to obtain a sustained spectrum 
oc E"^-^ from that energy up to Eankie [S]- 

In this paper we explore a different possibility recently proposed in [6]. Namely, the 
knee could be caused by new interactions with dark matter particles (x) in our galaxy. 
This hypothesis requires that the cross section a^^ becomes strong and breaks the cosmic 
protons of energy above the threshold Eknee in their way to the Earth. It is supported by 
the possible excess in the diffuse gamma-ray flux observed by MILAGRO |7] . Here we study 
other possible implications and also the conditions (on the cross section and the dark matter 
distribution) that have to be satisfled for this scenario to really work. 
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2 Other knees and the change in the spectral index 



In addition to protons, the cosmic ray flux contains atomic nuclei (He, C, Fe, etc.). The 
abundance of the different species at energies around 10^ GeV has been measured in extensive 
air-shower experiments. In particular, KASCADE [HI |9] has observed that these fluxes also 
seem to exhibit a knee at an energy that increases with the atomic number. It is obvious 
that if there is an energy threshold E^nee where the p-x cross section becomes strong and 
produces the proton knee, also a nucleus N will experience a similar effect. 

As a first approximation, if is negligible at energies below Eknee then the nucleus will 
start interacting when the energy per nucleon reaches that same threshold. This implies 
that the knee scales linearly with the atomic mass A: 

E^ee ^ A E,^ee • (2) 

One can take into account, however, that ap-^ is not zero aX E < E^nee (a step function), and 
that its raise with the energy may be better described by the power law 

ap^{E) CX E° {E< Eknee) • (3) 

In that case the atomic knee is moved towards the proton knee Eknee by two factors. First, 
the electric charge Z > 1 of a nucleus increases its interaction strength with the turbulent 
magnetic fields in our galaxy [TU] . This implies that the average time r (and distance L = cr) 
that it is travelling from the source to the Earth also increases 



r oc 



R-^ , (4) 



where R = E/Z is the rigidity and the index (3 at these energies may take values between 
0.3 (typical for a Kolmogorov spectrum of magnetic fluctuations) and 0.6 (as deduced from 
the study of fluxes of stable and unstable nuclei of energy below the TeV). The probability 
that a cosmic ray interacts along its trajectory grows with r. 



Pint ^ 1 - e^'^^x^x^ ~ % cr , (5) 

where = p^/m^ is the number density of dark matter particles (the average depth x of 
a trajectory is pL, and the nucleus interaction length is l/(o"Ar^n^)). Therefore, a larger r 
requires a smaller a^^ to obtain the same probability of interaction. The second factor is 
that the nucleus-dark matter cross section is larger than cip^. One can estimate that 

a (E)--l ^PX^E/A) Gp^ < 1 mbarn ; 
Nxy ) - \ ^2/3 ap^{E/A) ap^ > 1 mbarn . ^ ' 
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Therefore, the interaction probabihty at the nucleus knee will coincide with the one at the 
proton knee if 

E^^^^^Ek,,,,A^Z^ . (7) 

For large values of a we recover the linear scaling with the mass number, whereas taking 
a = 2 and (3 = 0.5, the helium, carbon, and iron knees would move from AEknee, ^^Eknee 
and 55Eknee to 2Eknee, 3i?fcneeand 5Eknee, respectiveljj^. Therefore, we conclude that the 
hypothesis under study here seems able to explain the knees for the different nuclei observed 
by KASKADE. 

Another important point to address is whether these interactions could also explain the 
constant spectral index —3.1 in the flux between the knee and the ankle. Let us suppose 
that with no dark matter perturbing its propagation the cosmic ray flux would have followed 
a dependence up to E ankle- This implies that the probabihty of interaction changing 

the spectral index from —2.7 to —3.1 must be 

/ jt; \ -0-4 

PintiE > Ek^,,) ^ 1 - (- . (8) 

We will show that this probability can be obtained if above i?fc„ee the proton-dark matter 
cross section keeps growing with a given exponent a': 

OC {E > Eknee) ■ (9) 

The probability of interaction for a cosmic ray along a trajectory of length / is 

p,„i(/,E) = l-e-'^-"'^'. (10) 

Now let us assume that / is a random variable with an exponential probability distribution 
w{l) and an average value (l) = L: 

w{l) = ie-'/^ . (11) 

This is justified if the cosmic ray propagation can be modelled by a diffusive process driven 
by irregularities in the galactic magnetic field [H]. The diffusion coeficient is D ~ oc 
and lengths larger than L will be suppressed exponentially. The probability of interaction is 
then 

P.ntiE) = [dl p.ntil, E) W{1) = 1 ^— — . (12) 

^ PX^X 

Since at energies above Eknee ^"px^x-^ grows fast larger than 1, a constant spectral index 
needs that 

/ ^ X -0.4 

(Tp^n^L OC . (13) 



E. 



knee 



^Intermediate values in the position of the knee would be obtained for larger values of a 
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If grows between the knee and the ankle hke E^''^'^", with /5 = 0.3-0.6, then we can 
explain a sustained spectral index of —3.1 in the cosmic ray flux. 



3 Constraints and predictions at neutrino telescopes 

The large probability of interaction required to produce the knee means that 

a,,n,L = ^ - 1 . (14) 

It is easy to see that this can only be achieved with very large values of the cross section 
and of the dark matter density p^. In particular, we know that as cosmic rays propagate 
from the source a significant fraction of them interact with interestellar matter. Since they 
collide with a hadronic cross section, the amount of matter that they cross must be a few 
g/cm^ (a hadronic interaction length is around 50 g/cm^ [1]). Therefore, on dimensional 
grounds x may have an impact on the cosmic ray propagation if ap^ grows above the mbarn 
and is larger than the density of interstellar gas. 

Let us first discuss the cross section. A possibility that has been extensively discussed 
during the past years is strong TeV gravity in models with extra dimensions [12]. If the 
fundamental scale is Md ^ 1 TeV, at larger center of mass energies the gravitational cross 
section will dominate over gauge boson exchange due to the spin 2 of the graviton. At 
lower energies (below the mass Mc of the first Kaluza-Klein mode) 4-dimensional gravity is 
unchanged, whereas between Mc and Newton's constant grows like a power law. An 
analysis of the proton-dark matter cross section at center if mass energies s > M|, in these 
models can be found in |13j. The scattering is dominated by eikonal processes where a 
parton carrying a fraction x of the proton momentum transfers a small fraction y of energy 
to X- The cross section at the parton level growj^ like gi+^Z"^ go it is larger and softer 
(involves larger distances and smaller energy transfer) for lower values of the number n of 
extra dimensions. This is easily understood because gravity dilutes faster with the distance 
for larger ra, if n < 2 it is a long distance interaction with a divergent total cross section. 

We find that if > 1 TeV only n = 1 seems to give the large values of ap^ required 
to explain the cosmic ray knee. Of course, we need in that case a mechanism to avoid 
bounds from supernova and from macroscopic gravity. For example, Giudice et al [I4j build 
a model where a warp factor adds a mass of ~ 10 MeV to the Kaluza-Klein excitations 
of the graviton without changing the basic properties of the 5-dimensional model at high 
energies. The effective model has a TeV fundamental scale and only one extra dimension up 



^Notice that growth of with the energy wiU also be determined by the parton distribution functions. 
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Figure 1: Cross section for an eikonal gravitational interaction in p-x collisions. 

to distances of ~ (10 MeV)^^. In Fig. 1 we plot the cross section in this setup. The average 
fraction of energy transferred hj a. E = 10® GeV proton to the dark matter particle in an 
interaction is (y) = 10~^. We have required a minimum of 1 GeV^, since the collision 
must break the incident proton. 

Let us now comment on the dark matter density that is required. Although it is thought 
that the total amount of dark matter in the galaxy is much larger than the amount of baryonic 
interstellar matter, the former would be distributed inside a spherical halo of ~ 200 kpc [15j, 
whereas the latter would be mainly within a disc of thickness ~ 6 kpc. As a consequence, 
the depth of dark matter crossed by a cosmic ray from the source (in the galactic disc) to 
the Earth would be similar or even smaller than that of baryonic gas (x/a/). For example, 
if we assume a constant ~ 0.3 g/cm^ (the local density near the solar system) then 

~ 0.1a;/A/- This value is clearly insufficient to explain the cosmic ray knee. 

We think, however, that one may consider scenarios where the density that cosmic 
rays face in their way to the Earth is substantially larger {x^ ~ 100 xjm)- In particular, 
the dark matter could be more clumped locally, in the regions where the cosmic rays are 
produced. It is also possible that cosmic rays spend a significant fraction of time inside local 
clouds of dark matter, trapped by stronger magnetic fields. They would also face a larger 
depth if the distribution of dark matter were more fiattened towards the galactic plane (a 
non-sperical halo) p^. This third possibility might be better accommodated if the galactic 
dark matter has two components (a heavier and a lighter one [Hj), as the proximity of the 
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Figure 2: Differential cross section da /dy (mbarn) in a v-p [or p-x) gravitational interaction, 
where y = {E^ — E^) / E^, is the fraction of energy lost by the incident neutrino (or transferred 
to x)- 

heavier component to the fundamental scale in these models would imply a larger x^X elastic 
cross section justifying an anomalous distribution. Notice that although the total amount 
of dark matter in the galaxy is well established, its distribution is underconstrained (the 
rotation curves can be fitted with different distributions [TBI [E] , and magnetic fields could 
also play a role PU]). 

We would like to discuss a last consequence that could be used to disregard the scenario 
under study in a model independent way. If gravity produces an increase in the dark matter 
interactions, it will also imply a large cross section for neutrinos above a certain energy 
threshold. More precisely, the center of mass energy when a neutrino of energy E^, hits a 
proton at rest is 

^^=^2^^^. (15) 

This coincides with the ^/s in the collision of a proton at the knee with a dark matter particle 
if 

= ^ Eknee ^ 10« GcV . (16) 

nip 

Therefore, in this framework we would not observe any neutrino events above El^ at z/- 
telescopes. Such a neutrino would have an interaction lenght of 1-10 m in rock, losing a 
y ^ IQ-^ fraction of its energy in each collision. In Fig. 2 we plot the distribution of y in 
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neutrino-proton interactions. 

For example, a cosmogenic neutrino [21] of 10^ GeV would interact with a mbarn cross 
section and lose around 1 TeV of energy {y 10~^) in each interaction with matter. Energy 
loss, however, is dominated by the (less frequent) collisions of larger y. We estimate that 
such neutrino could deposit a 10% of its energy every 1-10 km. This means that it would 
never reach a telescope from up-going directions. One would have, however, other possible 
signals there. The neutrino could reach the center of IceCube (2 km under the ice) [22] 
from large (inclined) zenith angles. It would start there a TeV hadronic shower every 1-10 
meters of ice, an event that would be clearly different from the muon bundle produced by 
an extensive air shower. 

Notice also that bounds on the v-p cross section from the non-observation of quasi- 
horizontal extended air showers [23] produced by cosmogenic neutrinos do not apply here. 
The reason is that in the interaction the neutrino will deposit just a small fraction y of 
energy, starting an air shower below the ~ 10*^ GeV threshold in these experiments. 

4 Summary and discussion 

Cosmic rays may be affected by new physics at the TeV scale. In particular, their interaction 
with dark matter seems interesting because of two generic reasons: ijthe dark matter particle 
is expected to be heavy, which provides larger center of mass energies {^2m^E) than in 
cosmic ray collisions with a nucleon, and ii) the standard interactions are expected to be 
weak, which makes the relative effect of new physics easier to detect. 

We have studied the possibility that the knee observed in the cosmic ray flux at energies 
around 10^ GeV may be caused by these interactions. In particular, models with extra 
dimensions and a TeV fundamental scale of gravity predict cross sections that grow very 
fast in the trans- Planckian regime. In p] it is shown that these interactions could break 
the protons and produce a diffuse flux of secondary gamma rays consistent with MILAGRO 
observations. Here we have found that the knee would also imply knees that scale with the 
atomic number for the other nuclei present in the cosmic ray flux. We have shown that the 
constant power law E~^'^ for the flux between the knee and the ankle could be correlated 
with a sustained growth oc E"' in the proton-dark matter cross section. 

We have also found, however, that the cross sections and the dark matter densities which 
are required are very large (mbarn and 100 times over the expected value, respectively). 
In particular, the size and the scaling of seem possible only in models with one extra 
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dimension. Such a model should be completed with a mechanism increasing the Kaluza- 
Klein masses like the one discussed in [H] to be consistent with observations. As for the 
dark matter, it may require a large local concentration near the sources of cosmic rays, 
or a flattened galactic distribution. We have speculated that this could be justified in a 
2-component model (see, for example, [T7]) where the heavier one may not be completely 
collisionless (the elastic scattering with another dark matter particle could have a center of 
mass energy just below the fundamental scale Md)- 

Notice also that the presence of a possible second knee in the cosmic ray spectrum at 
10^'^ GeV would not alter or invalidate the proposed mechanism. This second knee would 
have a different origin, the dominance in the spectrum of extragalactic protons [2l[[25]. The 
extragalactic comic ray flux would not be affected by the interactions producing the first 
knee because dark matter densities are smaller outside the galaxy and the trajectories at 
these energies are not a random walk [6] . 

In any case, the scenario that we have discussed has what we think is a model independent 
implication in neutrino physics that could prove it wrong: the absence of standard neutrino 
interactions above a threshold energy around 10^ GeV. Neutrinos above this energy would 
interact with protons with a mbarn cross section, losing a small amount of energy (around 1 
TeV) in each interaction. When a cosmogenic neutrino enters the atmosphere horizontally, 
the penetrating shower that it starts would contain an energy below the triggers in extended 
air shower experiments. However, in //-telescopes one could observe inclined events (down- 
going but from large zenith angles) where a very energetic neutrino starts a continuous of 
showers spaced by 1-10 m, a signal with no background within the standard model. The 
correlation between dark matter and neutrino physics is a generic feature in TeV gravity 
models. We find very appealing that a possible observational effect, the knee from cosmic 
ray-dark matter collisions, can be disproved (or supported) in the near future by another 
one, the observation of cosmogenic neutrino interactions at ICECUBE or ANITA |26j . 
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